Transport in the XX chain at zero temperature: 
Emergence of fiat magnetization profiles 



T. Antal, Z. Racz, A. Rakos 
Institute for Theoretical Physics, Eotvos University 
1117 Budapest, Pdzmdny setdny 1/a, Hungary 

G. M. Schiitz 

Institut fur Festkorperforschung, Forschungszentrum Jiilich 
52425 Jiilich, Germany 
(February 1, 2008) 

We study the connection between magnetization transport and magnetization profiles in zero- 
temperature XX chains. The time evolution of the transverse magnetization, m(x,t), is calculated 
using an inhomogeneous initial state that is the ground state at fixed magnetization but with m 
reversed from —mo for x < to mo for x > 0. In the long-time limit, the magnetization evolves 
into a scaling form m(x,t) = <&(x/t) and the profile develops a flat part (m = $ = 0) in the 
\ x /t\ < c(mo) region. The flat region shrinks to zero if mo — > 1/2 while it expands with the 
maximum velocity, Co = 1, for mo — > 0. The states emerging in the scaling limit are compared 
to those of a homogeneous system where the same magnetization current is driven by a bulk field, 
and we find that the expectation values of various quantities (energy, occupation number in the 
fermionic representation) agree in the two systems. 
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I. INTRODUCTION 



Transport in intcgrable systems is anomalous. This 
means that a harmonic lattice (or a transverse Ising 
chain) cannot support an internal thermal gradient when 
the two ends of the system are kept at different temper- 
atures [fUH]. The temperature profile that forms is flat 
everywhere apart from the neighborhood of the bound- 
aries, and the energy flux is clearly not proportional to 
the temperature gradient inside the sample. The origin 
of this anomaly is the fact that currents in integrable 
systems (energy currents in the above examples) often 
emerge as integrals of motion thus causing the transport 
coefficients to be divergent or ill defined J3j. 

The existence of flat density (temperature) profiles in 
the presence of (energy) currents points to the important 
role played by the homogeneous, current-carrying states 
in intcgrable systems. In previous papers ||,||, we sug- 
gested that a method of constructing such states at zero 
temperature is to add to the Hamiltonian an appropri- 
ate current with a Lagrange multiplier and then find the 
ground state of the system. This method can be applied 
with relative ease to models such as the transverse Ising 
and isotropic XY (so called XX) models, and one can 
calculate the effect of energy current on experimentally 
measurable correlations. 

Since there are physical systems that are rather well 
approximated as XX models Hj, and since the effect of 
currents should be measurable in inelastic neutron scat- 
tering experiments j7|, one should carefully examine the 
validity of the assumptions underlying the calculations 
of the correlations. The Lagrange multiplier method is 
based on the assumption that the homogeneous, current- 



carrying states are the same whether they were induced 
by boundary conditions or by bulk driving fields. This 
is not an obvious assumption and its validity should be 
checked by first solving the boundary condition problem 
and then comparing the results with those obtained from 
the Lagrange-multiplier method. Unfortunately, the ex- 
ecution of this program runs into difficulties since the 
boundary condition problem is nontrivial and we cannot 
provide an exact solution for any of the interesting cases. 

One can make some progress, however, in a closely 
related problem. Namely, current-carrying states gener- 
ated by inhomogeneous initial conditions can be studied 
and one can ask whether those states were describable by 
the Lagrange multiplier method. In order to illuminate 
the new problem, let us imagine the original problem as 
finding the steady state of a brick with two of its opposite 
sides kept at different temperatures. In the new problem, 
we have two semi-infinite bricks at different temperatures 
and, at time zero, the bricks are joined at the ends. If 
the bricks were integrable systems then one would ex- 
pect that the time-evolution would produce an interest- 
ing temperature profile with a flat part around the joining 
point and it would be meaningful to compare the pro- 
file as well as various local expectation values to those 
obtained from the Lagrange-multiplier method. This is 
what we shall do in this paper for the XX chain after 
making the following simplifications. First, we note that 
the transverse magnetization, to, is a conserved quantity 
in the XX model. Since this is a much simpler quantity 
to treat, we shall study the transport of to instead of the 
transport of the energy. Second, we restrict our study 
to the zero-temperature properties of the system. Third, 
we shall use a simple inhomogeneous initial state pasted 
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together from two zero-temperature homogeneous pieces 
with m(x < 0, 0) = mo and m[x > 0, 0) = —mo. 

As a result of the above simplifications, the problem of 
time-evolution becomes solvable and the formation of a 
flat magnetization profile can be observed. We find that 
both the magnetization profile and the fermionic two- 
point correlations agree with the results obtained from 
the Lagrange multiplier method, thus indicating that the 
method does work for obtaining the homogeneous current 
carrying states. 

We start (Sec|nj) by defining the system and construct- 
ing the initial state in terms of fermionic operators. Next, 
the time evolution of the magnetization is calculated and 
the scaling form, m(x, t) = $(x/t), of the long-time limit 
is discussed (Sec. III). Finally, we show that the local 
properties of the scaling limit can be described using the 
homogeneous, current carrying states obtained from the 
Lagrange multiplier method (Sec{n~ 



transform provides us with the time dependence of a lo- 
cal Fermi operator, c„, as a sum of Bessel functions of 
the first kind 101 



(4) 



The above expression is the starting point for our calcu- 
lations since we wish study the local transverse magneti- 
zation, ,Sj l5 that can be expressed through the local Fermi 
operators as 
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Our aim is to find the time evolution of the expectation 
value of st 



m(n,t) = (<p\s*(t)\<p) 



(6) 



II. THE MODEL AND THE CONSTRUCTION OF 
THE INITIAL STATE 



using an initial state \<p) that is expected to produce a 
magnetization current. The simplest such state has a 
symmetric step in the initial magnetization 



The system we investigate is the XX model defined by 
the Hamiltonian 



H' 



Y, h n 



J > ( s n s n+l + s n s 



y s v ) 



(1) 



where the spins, (a = x, y, z), are 1/2 times the Pauli 
matrices situated at the sites of a one-dimensional peri- 
odic lattice [n = 0, ±1, ±(N - 1),N ; s% 



N+l 



-N+ll- 



The coupling, J, is set to J = 1 throughout this paper 
and N — > oo is assumed in most of our calculations. 
The z component of the total magnetization, M z = 
s n' ^ s conserved in this system and, using the conti- 
nuity equation for M z , a local magnetization current can 
be defined lb|: 



J 



M 



'n*n+l 



(2) 



It can be easily shown that the global magnetization cur- 
rent, J M = ^2 n jn , commutes with the Hamiltonian (|l|) 
and thus one can classify the eigenstates of .ff xx by the 
global current flowing through the system. Below we 
shall study how the initial conditions can generate states 
that are locally homogeneous and carry a finite magneti- 
zation current. 

The time-evolution that follows from a given initial 
condition can be calculated for the XX model since, us- 
ing a Jordan- Wigner transformation, this model can be 
transformed into a set of free fermions 101 



H xx = _J2 cosk 6 t 6fc 



(3) 



where bk is the annihilation operator of a fermion of wave- 
number k. It follows that the time evolution of the bk op- 
erators is known bk(t) = bk(0) exp(ztcos k) and a Fourier 



/ l zfnM \ / m for —N < n <0 , . 

= { -m for l<n<N. (7) 

The above condition, however, does not specify the initial 
state uniquely. In order to define \ip) precisely, we shall 
construct the state of the 2N spins as a direct product 
of states of two systems of N spins. The left system con- 
tains the n < spins while the right one is built from the 
rest (n > 1). In order to get close to a zero-temperature, 
current-carrying state we choose the lowest energy state 
under the above conditions (Q) . Namely, both half-chains 
are chosen to be in their ground states with the magne- 
tizations being mo and —mo in the left and right, respec- 
tively. Another reason for choosing such a state comes 
from trying to avoid possible energy currents in the sys- 
tem. In the N — > oo limit, the above state is a ground 
state of the model at fixed m = ±mo and the two homo- 
geneous sides have the same energy. Thus one expects 
that only magnetization current will flow in the course of 
the time-evolution of the system. 

A formal definition of \ip) with the above properties 
can be given as follows. Let us define annihilation oper- 
ators on the left and on right by the Fourier transforms 
of the fermionic operators, c„, on the appropriate sides 
of the system 
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As can be shown, these operators are also fermionic op- 
erators with the anticommutation relations: {Lk, L+} = 
{Rk, Rq} = $kq, and the initial state, \<p), can be written 
as 
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fe=-fc_ 



fc'=-fej 



where | J, ) is the state with all spins down and k± = 
7r(i dz mo). In the following m > will be assumed 
without restricting the generality of the arguments. 



III. NUMERICAL AND ANALYTICAL 
EVALUATION OF THE TRANSVERSE 
MAGNETIZATION 

Once the initial state \ip) is given, the time evolution 
of s z n is expressed through initial correlations 
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It follows now from the special construction of \<p) that 
the terms with j and I on different side of the origin van- 
ish. Furthermore, the terms where I — j is odd cancel 
since the expectation value is real. The calculation of 
((p\cf (0)cj(0)\ip) is carried out by using the explicit form 
(eqs. p|) for \<p) and the result is as follows 
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j even ^ 
j = • 
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The upper (lower) sign is valid for both s and j being in 
the left (right) part of the chain. Since the above expres- 
sion depends only on the even integer s — j, the nota- 
tion simplifies considerably by introducing I = (s — j)/2. 
Treating the I = term separately and using an iden- 
tity [Q for the Bessel functions, J^n^n = 1; equations 
( |l0| ) and ([ll|) yield the following magnetization profile 
for n > 1: 



m(n, t) 



n-l 
j = l-n 



n-l 
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(12) 
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The n < values can also be calculated and they show 
the symmetry of the initial state, m(n,t) = —m[—n + 
l,t). The above expression for m(n, t) is the basis for 
both numerical and analytical results presented below. 



A. Numerical results 

The infinite sum ( p"2| ) defining m(n,t) converges well 
since the main contribution comes from terms with small 



I. The numerical evaluation of m(n,t) poses a problem 
only for mo <C 1/2 when the scaling limit, n — > oo, 
t — > do, with n/t r*> 1, can not be readily reached due 
to the increasing number of terms that must be taken 
into account when evaluating the sums in (p"2|). 

The results from the numerical analysis (see Fig.|l|) 
show that the magnetization profile evolves into a scaling 
form in the large-time limit 



m(n,t) w $(n/i) 



(13) 



with the following scaling function (analytically obtained 
in Sec, ill B ) giving an excellent fit to the data: 



$(u) = = 
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- arccos(f ) 



for < « < cos(7rmo) 

for cos(7rmo) < v < 1 (14) 

for 1 < v 
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FIG. 1. Magnetization profile for initial conditions 
mo = 0.5 and mo = 0.2. The large-time limits approach 
the scaling curves (solid lines) given by equation (0). 

As one can see, the magnetization profile does de- 
velop a flat part that expands with a finite velocity 
c = cos7rmo. The velocity c(mo) decreases with increas- 
ing mo and it diminishes for mo —►1/2. 

The conservation of magnetization implies that the 
magnetization current has a scaling form, j M (n,t) w 
^>(n/t), as well, and the scaling function can be obtained 
from ( |l4| ) using the continuity equation together with 
the boundary condition j AI (n — > ±oo,t) = 0. We find 
the scaling function 




for < v < cos(7rmo) 
for cos(7rmo) < v < 1 
for 1 < v , 



(15) 



shown in Fig.p 
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FIG. 2. Scaling functions for the magnetization-current for 
initial conditions mo = 0.5 and mo = 0.2 as given by equation 
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B. Analytical results 



The case v < 1 is complicated by the fact that the 
limit n — > oo does not exist for ^' n [v) due to the increas- 
ing frequency of oscillations in ( jig ) . It should be noted, 
however, that we are interested in the integral of $>'(v) 

§(«) = lim $ n (u) = lim / , (20) 

n — >oc n — >oo Jq 

and the n — > oo limit of the above integral does exists. 
The value of the integral can be obtained by replacing the 
fast oscillating cos 2 term [see (|l9|)] by its average value 
1/2. One then obtains 

$>(v) = arcsinu , (21) 

in agreement with ( p"4| ) for m = 1/2. 

The evaluation of $(u) for to 7^ 1/2 is more involved 
but, essentially, it follows the above steps and yields the 
general expression (|L1| ) . An outline of the calculation and 
some intermediate results are given in Appendix A. 



Let us begin by deriving the analytical guess obtained 
from numerics (|l4|) by treating the case too = 1/2. This 
limiting case is especially simple since sin (27ttooZ) = 
and hence one is left with the following expression: 

™(M = 4 ]T 4w = -J + E J 'w- ( 16 ) 



j=l-n 



j>n 



We shall evaluate the scaling limit of the above sum by 
calculating the integral of the large n limit of the discrete 
derivative of m(n,t) defined as |ll|] 

& n (v)=t [m(n+l,t)-m(n,t)] n/t=v = •%(-) ■ 

(17) 

The cases v > 1 and v < 1 must be treated separately 
since the asymptotic properties of the Bessel functions 
change at v = 1 [Eof. Namely, one finds for v > 1 
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(18) 



while for v < 1 
<(«) ~ —± 

7TVl ^ 



( \/ 1 — v 2 — v arccos v) 



(19) 



As one can see, the derivative approaches zero exponen- 
tially for v > 1 (outside of the 'light cone') and thus $(u) 
is constant in this region. The constant is determined by 
the boundary condition $(00) = —1/2 and thus we have 
$(V) = -1/2 for v > 1. 



IV. COMPARISON WITH THE LAGRANGE 
MULTIPLIER METHOD 

A consequence of the emergence of scaling, m(n,t) w 
$(n/f), is that the magnetization becomes locally flat, 
m(n + l,i) — m(n,t) ~ — > 0, in the scaling limit. 
Since these locally homogeneous segments of the system 
become arbitrarily large for t — ► 00, one would like to 
ask if they could be described as a homogeneous, current- 
carrying state of an infinite chain as obtained from the 
Lagrange multiplier method || . 

In order to formulate the question more precisely, let 
us note that, for a given n/t = v and in the neighborhood 
of n, the scaling state of the system may be characterized 
by the local values of the conserved quantities and their 
fluxes, where the expectation value of a local operator 
a n (t) is defined in the scaling limit as 



= lim (ip\a n (t)\ip) 

n.t— >oo 



(22) 



n/t—v 



The simplest conserved quantities are the magnetization, 
to = (s z ) v = $(?;), and the energy density, e = (h xx ) v 
with the corresponding fluxes being the magnetization 
current, j M = (j M ) V i and the internal energy flux, 
j 1 = (j^v, where 



(23) 



The energy density is homogeneous in the initial state 
(apart from the local perturbation of the domain wall at 
the origin) and it can be shown that this homogeneity re- 
mains intact in the long-time limit. Thus e is determined 
by the initial condition and, in the scaling limit, can be 
expressed through to and j M 
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e = e(m,j ) = cos(7rmo) = 

cos [7ijm| + arcsin(7rj M )] . (24) 

Note that the homogeneity of the energy density implies 
a vanishing energy flux, j 1 — 0, for all values of v and 
so equation (|24|) can be viewed as an equation of state 
connecting the relevant densities and fluxes. 

It is clear that (m,j M , £,i J ) gives only a partial char- 
acterization of the state since there are infinitely many 
conservation laws in the XX model. Nevertheless, one 
may assume that (m, j M , e, j 1 ) were sufficient to describe 
the states emerging in the scaling limit. Then one can ask 
if the properties of these states were the same as those 
of the homogeneous, current-carrying state of an infinite 
system where the values of m,j M were fixed by in- 
troducing conjugate fields (Lagrange multipliers) and by 
finding the ground state, \ip), of the following Hamilto- 
nian 

H = H XX - hM z - X M J M -XiJ 1 . (25) 

The reason for this identification is that the scaling 
emerges as an evolution from an initial state of mini- 
mal energy at a given m — mo- We assume here that 
local states of the scaling limit remain zero temperature 
states in the sense that they have minimal energy density 
at the given values of m,j M , and j 1 . 

Since m, j M , j 1 are set by adjusting the Lagrange mul- 
tipliers, a necessary condition for \ip) having the local 
properties of the scaling states is the equality of its en- 
ergy density as a function of m, j , j 1 to the energy den- 
sity e given by (^4|) . Verification of this equality can be 
viewed as checking whether the equation of states of the 
two systems were identical. 

In order to obtain (iplh^ , we use the results from 
our previous work J5j where the problem of energy cur- 
rent in the transverse XX model was investigated. In 
that work, we studied the ground-state properties of the 
following Hamiltonian 

U = H XX - hM z - X E J e ■ (26) 

with the total energy flux, J E , being a sum of the fluxes 
of internal and of magnetic energy 

J E = J 1 - h,J M . (27) 

We found that the j E , h phase diagram of this model has 
a region (called region (g) in Fig. I of ||, see Fig.|| of the 
present paper) where the ground-state expectation value 
of the flux of internal energy is zero ((iplJ^ip) = 0) and 
so the j 1 — condition is automatically satisfied. The 
magnetization flux in region @ has been calculated in 
||. For the case of hXE < that is needed in the setup 
we use (mo > 0), it has the following form 



and the calculation of the densities of magnetization and 
of energy is also a simple exercise with the following re- 
sults 

m = — (arcsin/i — arcsinA^, 1 ) — isign(fe) , (29) 



Whtrw) = --V^k 2 ■ (so) 

Eliminating h and Xe from (|2|j2^ and ^) produces then 
(p4|) , i.e. (^>|/i* x |^) = e(m,j M ). This result indicates 
that the local states derived as a scaling limit can indeed 
be obtained with the help of the Lagrange method. 
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FIG. 3. Magnetic-field (h) - energy-flux (j E = j 1 - hj M ) 
phase diagram of an XX model where the energy flux is 
driven homogeneously in the bulk [5]. The dashed lines in 
region © correspond to states of constant energy with the lo- 
cation of the dashed lines determined by the initial conditions. 
The two filled squares represent the m = ±mo regions while 
the filled circles describe identical states corresponding to the 
plato (m=0) in Fig.l. Moving along the dashed lines provides 
the local states (m,j ,j = 0) in the transient region. 

Further confirmation of the validity of this sugges- 
tion comes from the studies of the scaling limit of the 
fermionic number operator (rik)v in momentum space, 
which can be defined as follows: 

(n k ) v = J2e- ikm Hm {ip\c+ t+m {t)c vt {t)\ip) (31) 

ra 

This quantity can be calculated for the special case of 
too = 1/2 (see Appendix B) and one finds that 

(nk)v = 
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J 1 if — arccos('y) + tt/2 < k < arccos(v) + tt/2, , ■. 
\0 else. {62) 

Since nk is a projection operator with eigenvalues 0, 1 on 
the eigenstates of H xx this result shows that the system 
evolves locally (in the scaling sense described by (gg)) 
into an eigenstate of the XX chain. The occupied levels 
are the same as those obtained by the Lagrange multiplier 
method for the given set of (m, j M , j 1 = 0). This can be 
interpreted as the local scaling states being equivalent 
to the ground states found in the Lagrange multiplier 
method. 

In order to demonstrate the value of the equivalence 
suggested by the above considerations, let us show how 
simple is the derivation of the magnetization profile if we 
assume that the Lagrange multiplier method can be used 
to describe the scaling limit. 

The calculation is based on the observation that the lo- 
cally flat magnetization [m(n +l,t) — m(n, t) ~ t -1 <fr' — * 
0] implies that one can take the continuum limit of the 
equation describing the conservation of magnetization 



d t m + d xJ M = 



(33) 



Next, using the fact that the energy density is constant 
in both space and time one finds the states that can be 
used to describe the scaling limit for a given initial mag- 
netization. It follows from eq. (|30|) that these states are 
located along h — constant lines in the j E , h phase dia- 
gram as shown in Fig.|§|. Along this path, the current j 
depends only on the local magnetization 
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j (x, t) = F(m) — — sin [n(m — \m\ 

7T 



(34) 



where F(m) is obtained from (g8|) and ( p9| ) using the 
fact that h is given by the initial energy density through 
(|30|). It follows then that the continuity equation (|3^ ) 
has a scaling solution [including the solution m = 0] and 



that the magnetization profile, m(x,t) = &(x/t) can be 
obtained from the following equation 



dF 
dm 



x 
1 



(35) 



Taking into account the limiting behavior, <$>(z — > 
±oo) = ±m , the straightforward steps described above 
lead to the results ([l4]) obtained from a much more com- 
plicated calculation treating the dynamics of the system. 



V. FINAL REMARKS 

We have shown that, in the zero-temperature XX 
model, an initial state with a step-like magnetization in- 
homogeneity evolves into a scaling state that is remark- 
able for its magnetization profile (emergence of a flat 
part) and for the existence of locally homogeneous, cur- 
rent carrying states. We have also provided evidence that 
the local states of the scaling limit can be obtained by 
the Lagrange multiplier method. The possibility of ap- 
plying the Lagrange method would simplify the calcula- 
tion of transport significantly, thus it would be important 
to test the generality of our results. The simplest ques- 
tions to ask are whether the results could be extended 
to finite temperatures and whether the calculation could 
be generalized to the transport of energy as well. It ap- 
pears that these questions can be answered and work is 
in progress along these lines. 



ACKNOWLEDGMENTS 

This work has been partially supported by the by the 
Hungarian National Science Foundation (Grants OTKA 
T-019451 and T-017493). 



APPENDIX A: MAGNETIZATION PROFILE 



We are considering the case mo 7^ 1/2 and study the scaling limit of the discrete derivative & n {v) defined in (|17|): 
*n(«)= 1 + M) - m{n, t)] n/t=v 

2_rnon j2 g sin(Wr) + 



(Al) 



ijt—v 



In order to take the scaling limit, we divide the sum into three parts where the indexes of the Bessel functions are 
positive. Using the limiting properties of the Bessel functions, one obtains the following form (n — > 00 and < v < 1): 



TT 2 y/l — V 2 
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(n tt 
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A* = lim 



sin(2mo7rr) 



n^oo 27rr 

sin(2mo7rr) 

r C< 

27rr 



cos 



2rv\ 7r 



Tl«± — ) 



7Ta/1 — (f ± 2rv/n) 2 



n 7r 

— 7(f) 2rarccost> 



7T-S/1 



(A3) 



where 7(2;) = Vl — x 2 — x arccos x. We can drop the second sum since it vanishes with n — > 00. Using trigonometrical 
identities, the calculation is reduced to evaluating the following sums (for < x < 2tt) 



E 



sin(rcc) n — x 



and E 



cos(rcc) 



--ln(2(l - cosx)). 



(A4) 



and one finds 



, . . JO for v < cos mo7r, 

*«->oo(. u ) - |_ cos 2 (~7(t>) - 7r/4) for w > cosm 7r. 



(A5) 



Just as in Sec. Ill B one can see that <£>' ^= lim, woo since the differentiation and the scaling limit cannot be 
exchanged. However, one can obtain $ from ( |A5| ) by simply integrating it with the substitution cos 2 — * 1/2 as it is a 
rapidly oscillating function of v in the scaling limit 



<&(v) = lim / & n 



(y)dy 



for v < cosmo7r, 

rv dy _ _ m _|_ i arccosw f or v > C0Smn7T. 



(A6) 



APPENDIX B: FERMION OCCUPATION NUMBER IN THE SCALING LIMIT 

In a translationally invariant system the expectation value of the fermionic number operator (rik) = ( b\bk ) in 
momentum space is expressed in terms of the real-space creation and annihilation operators as the Fourier transform 
of the two-point correlation function 

K>=E e " fem <4+™9> (bi) 

m 

where the site number j is arbitrary (because of translational invariance). In order to investigate this expectation 
value in the scaling regime we choose j = vt and take the stationary limit t — ► 00 in the expectation value before 
performing the summation over m. In this way we probe the scaling state of the system in the region of the shifted 
origin located in the point vt where the state becomes translationally invariant in the long-time limit. 
At time t the correlation function (cj +Tn Cj ) for the initial state with mo = 1/2 is obtained from (||) as 

00 

(4+™ c ,> = H) m E J s+™+;(*)-W*) 

(—i) m t 

= (— i) Jm+jJj H Yj {Jm+jJj+l - Jm+j+lJj) ■ (B2) 

Using the asymptotic expression of the Bessel functions Eo) 
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Jvt-\-m if) 



TTtVT 



= cos (ty/ 1 — v 2 — (vt + m)a — — ) 



(B3) 



with a — arccos v and utilizing the addition theorems for trigonometric functions we thus obtain 



lim ( c. 



t 

vt-\-m 



(t)c vt (t)] 



irm 
1 

2^ 



a+n/2 



dp e 



iprn 



-tt/2 



The Fourier transformation finally yields 



and hence (|32|). 



a+n/2 



./2 



a+n/2 



dp S(p — k) 



-a+Tr/2 



(B4) 



(B5) 
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